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Abstract 


Drought hazard is one of the natural features of the planet that can occur in all climatic zones so 
that in some cases it causes irreparable damage. Due to the importance of drought, this study 
aimed to analyze the time series of agricultural drought in Darab city over 20 years (2000-2020) 
using remote sensing and satellite images in the Google Earth engine system. Satellite imagery 
including 460 images of land surface temperature (LST) and vegetation (NDVI) MODIS of Terra 
satellite was used to calculate the temperature condition index (TCI) and vegetation condition 
index (VCI). Moreover, using the recorded rainfall data of the Darab synoptic station, the SPI 
index was calculated by MATLAB software in different time intervals. The results show that 
according to the extreme dryness class and based on the TCI index, the 2000 with an area of 
225.46 square kilometers and also according to the VCI index, the 2013 year with an area of 
280.80 square kilometers, had the most dried land area in Darab city. As a result, comparing the 
numerical value of SPI index for each of the years of the period under study with the amount of 
non-drought areas obtained from TCI and VCI indices obtained from satellite images, the highest 
correlation coefficient of 0.76 was observed between the 12-month SPI and the VCI index, which 
indicates the VCI satellite index as the optimal index indicating the drought situation in Darab 
city. Another outstanding result of this study is that the use of remote sensing data and Google 
Earth Engine System to monitor and investigate drought in areas that do not have observational 
data from land surveys is very useful. 
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حكيده 


مخاطره خشکسالی ازجمله ویژگی‌های طبيعى كره زمين محسوب مىشود که امكان رخداد آن در تمام 
مناطق آب و هوايى وجود دارد؛ به‌طوری‌که در پاره‌ای موارد خسارات جبران‌ناپذیری را ايجاد می‌کند. با 
توجه به اهميت خشکسالی» اين پژوهش با هدف تحليل سرى زمانى خشکسالی کشاورزی شهرستان 
داراب طى یک دوره ۲۰ ساله (۱۳۹۹ - ۱۳۷۹) با استفاده از سنجش ازدور و تصاوير ماهواره‌ای در محيط 
سامانه گوگل ارث انجین انجام شد. تصاویر ماهواره‌ای شامل ۶1۰ تصوير از محصولات دمای سطح زمین 
LST)‏ و پوشش گیاهی (NDVD‏ سنجنده مودیس ماهواره ترا می‌باشند که از آن‌ها جهت محاسبه 
شاخص وضعیت دمایی TCI)‏ و شاخص وضعیت پوشش گیاهی (VCD‏ استفاده شد. همجنين با استفاده 
از داده‌های بارش ثبت شده ایستگاه سینوپتیک داراب شاخص SPI‏ به‌وسیله نرم‌افزار MATLAB‏ در 
بازه‌های زمانی مختلف محاسبه شد. نتایج به دست آمده حاکی از اين است که مطابق با طبقه خشکی بسیار 
شدید و بر اساس شاخص TCI‏ سال ۱۳۷۹ با مساحت ۲۲۵/۶۱ کیلومترمربع وسعت و همچنین نيز بر 


اساس شاخص VCI‏ سال ۲ با مساحت ۲۸۰/۸۰ کیلومترمربع وسعت. بیشترین مساحت خشکی را در 
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۱۳ جغرافیا و مخاطرات محيطى شمارة چهلم 


شهرستان داراب داشته‌اند. درنتیجه مقایسه بين مقدار عددی شاخحص SPI‏ برای هر یک از سال‌های دوره 
موردبررسی با میزان مساحت‌های فاقد خشکسالی حاصل از دو شاخص TCI‏ و VCI‏ به دست آمده از 
تصاویر ماهواره‌ای. بیشترین مقدار ضريب همبستگی به میزان ۰۸۷۲ بين SPI‏ دوازده‌ماهه و شاخص VCI‏ 
مشاهده شد که اين مطلب نمایانگر شاخص ماهواره‌ای: VCI‏ به‌عنوان شاخص بهينه نشان‌دهنده وضعیت 
خشکسالی در شهرستان داراب است. 


۱ - مقدمه 


خحشکس الی یکی از مخرب‌ترین خطرات طبیعی کره زمین است که میلیون‌ها فرد را در زمینه‌های مختلفی از قبيل 
امنیت غذايى» مهاجرت و ضررهای اقتصادی تحت تأثیر قرار می‌دهد (کولکارنی و همکاران e‏ ۰ و به‌صورت 
مداوم از چندین ماه تا چندین سال اتفاق می‌افتد كه در طی آن میزان آب مورد نياز مناطق کاهش می‌یابد (ویاس و 
همکاران أ ۵ به‌طوری که به‌عنوان جدی‌ترین مشکل در جهت تبعات اجتماعی محسوب می‌شود (هزایمه و 
حسن ‏ ۲۰۱۱) بااین‌حال هیچ تعریف پذیرفته شده و یکسان جهانی برای خش كس الى وجود ندارد (توولدبیرهان 
سیگ ۲۰۱۹) و ازانجایی‌که مشاهده تأثیرات مخرب آن به زمان طولانی نياز دارده نحوه مقابله با آن نیز کمتر 
موردتوحه واقع شده | (عليقلىنيا» (AVY AYIA‏ واضح است که خشکسالی دارای روندی آهسته و در حال 
تکامل است که با کمبود میزان بارش شروع شده و سپس منجر به از دست رفتن رطوبت خاک و افزايش دمای سطح 
کشاورزی مناطق دارد (شن و همکاران ۹ خشکسالی کشاورزی زمانی به وقوع می‌پیوندد که تأثیرات ناشی از 
خشكسالى بر محصولات کشاورزی قابل روئیت باشد (نعمتى و همكاران» ٤ AYIA‏ كه در این زمينه داده‌های 
حاصل از تصاویر ماهواره‌های سنجش ازدور قادرند به‌طور مداوم توزیع مکانی و زمانی از رشد محصولات و وضعیت 
خحشکسالی را ثبت نمايند (ژو و همکاران" (Yeye‏ بنابراین حشکسالی را می‌توان بر اساس تغییرات مشاهده شده در 
سلامت پوشش گیاهی و همجنين پوشش اراضى سطح زمين با استفاده از داده‌های سنجش‌ازدور مورد ارزيابى قرار 
داد SU)‏ و چودور“ A IAV‏ نعمانی و همکاران“ ۹ سیلئوس و همکاران أ (Yer‏ بدين جهت سنجثازدور 
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سال دهم بررسى سرى زمانى خشكسالى كشاورزى شهرستان داراب .... WV‏ 


به‌ویژه با كسترش شاخصهاى هدفمند wl‏ سابقه‌ای طولانى جهت سنجش و نظارت بر ویژگی‌های خشكسالى و 
تأثيرات bol‏ بر روى زمين دارد (آقاکوچک و همکاران؛ ۲۰۱۵). وانگ" و همكاران (۲۰۱۸) با استفاده از معادلات 
جند جمله‌ای» شاخص نرمال شده تفاضل پوشش كياهى (NDVI)‏ و دماى سطح زمين» شاخص ETVDI‏ را جهت 
بررسی خش كس الى مناطق ابداع و معرفى نمودند. جی* و همكاران (۲۰۱۸) با استفاده از داده‌های سنجش ازدور و 
جندين شاخحص خشکسالی. وضعيت خشكسالى سيجوان جين را موردبررسى قرار دادند. نتايج مطالعه آن‌ها حاكى از 
رخ داد خشکسالی شديد بين سال‌های ۲۰۰۰ تا ۲۰۱۱ در اين ناحيه مىباشد. کیو" و همكاران (۲۰۱۹) با استفاده از 
داده‌های ماهواره‌های سنجش زدور و شاخص‌های مختلف خش كسالى به تهيه نقشه نوع كشت مطابق با وضعيت 
خش كس الى برای كشورهاى واقع در شاخ افریقا اقدام نمودند. فنگ" و همکاران (۲۰۱۹) به‌وسیله داده‌های 
سنجش‌ازدور موديس و TRMM‏ و همجنين نيز استفاده از جندين شاخص محاسبه خشکسالی, نقشه‌های خشكسالى 
ناحيه جنوب شرقى استراليا را در طى فصول مختلف تهيه نمودند. فيصل” و همكاران (۲۰۲۰) به بررسى خشكسالى 
منطقه يايو غربى در اندوزى با استفاده از تصاوير سنجش زدور پرداختند كه نتيجه مطالعه آنها حاكى از قابليت 
بالاى سنجش زدور در بررسی وضعيت JUSS‏ می‌باشد. بيك '' و همكاران (۲۰۳۰) به ارزیابی خشكسالى 
هواشناسی و حشکسالی كشاورزى در حوضه رودخانه چیترال "" كابل با استفاده از داده‌های سنجش‌ازدور و جندين 
شاخص خشکسالی اقدام نمودند. نتيجه كلى مطالعه آنها بیانگر اين است که خش كسالى را مىتوان با ادغام منابع 
مختلف داده به‌درستی ارزيابى نمود. توولدبیرهان سيك '! و همكاران (۲۰۲۰) به مقايسه شاخصهاى مختلف 
خشکسالی هواشناسى و خشكسالى كشاورزى در اتيويى پرداختند. آن‌ها بیان نمودند كه جهت برنامهريزى بهتر در 
مقابله با خشکسالی» بايد تفاوتهاى بين خشکسالی كشاورزى و هواشناسی هر منطقه مشخص شود. ژو " و 
همکاران (۲۰۲۰) با استفاده از تصاوير ماهواره‌ای Sentinel-2‏ و مودیس و همچنین نيز استفاده از شاخص‌های 
NDVI‏ و VICI‏ نقشه وضعیت خشکسالی دشت گوانژونگ در جين را تهیه نمودند. سلطانی‌فر و همکاران (۱۳۹۸) 
در مطالعه‌ای شهرستان کرمانشاه را با استفاده از تصاویر لندست به لحاظ حشکسالی مورد پایش قرار دادند. نتایج 
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مطالعه آن‌ها نشان داد که تصاوير لندست داراى قابليت مورد نياز جهت پایش خش كسالى می‌باشند. زرگران و 
همكاران (۱۳۹۸) به يايش خش كسالى روزانه استان كرمانشاه با استفاده از داده‌های پوشش كياهى حاصل از شاخص 
NDVI‏ سنجنده موديس و داده‌های بارش شبكه INTERIM‏ يرداختند. وجود همبستكى قوی و معنىدار بين تغييرات 
بارندگی داده‌های INTERIM‏ و شاخص NDVI‏ از نتايج مطالعه آن‌ها می‌باشد. سليمانى ساردو و زارعى (۱۳۹۸) با 
استفاده از توليدات سنجنده موديس و شاخص بارش استاندارد شده SPI‏ حوضه آبخیز جازموريان را به لحاظ 
حش كسالى مورد پایش قرار دادند. نتایج مطالعه آن‌ها بيانكر اين است که بخش اعظمى از محدوده موردمطالعه در 
كستره كلاسهاى خشكسالى خيلى شديد و شديد قرا دارد. با توجه به موارد ذكر شده و ازآنجايىكه شهرستان داراب 
همجون بسيارى از مناطق كشور طى سالهاى اخير با مشكل مخاطره خشكسالى مواجه بوده و از این حيث خسارات 
بسیار زیادی در شهرستان به‌ویژه در بخش کشاورزی به وجود آمده و معضلاتی را ایجاد نموده است. لذا هدف از 
انجام اين مطالعه پایش و بررسی سری زمانی بروز خشکسالی شهرستان داراب با استفاده از تصاوير ماهواره‌ای و 
محصولات سنجنده مودیس ماهواره ترا" در محيط برنامه‌نویسی سامانه تحت وب گوگل ارث انجین " می‌باشد تا 
بدین طریق ضمن شناسایی نواحی تحت SU‏ خشسکسالی در شهرستانء به‌عنوان راهنمایی جهت برنامه‌ریزی 
دست‌اندرکاران امر در زمینه مقابله با مخاطره خشکسالی و پیامدهای آن مورداستفاده قرار گیرد. سامانه گوگل ارث 
انجین در طی سال‌های اخیر فرایندهای پردازش تصاوير ماهواره‌ای را سهل نموده و بدين لحاظ مورد استقبال بسیاری 
از پژوهشگران قرار گرفته اسست. از قابلیت‌های گوگل ارث انجین می‌توان به کاهش بسيار زياد زمان فرایندهای 
پردازش تصاوین پوشش وسیع و امکان دسترسی به تصاویر و داده‌های رایگان کل کره زمین» امکان ایجاد سرى 
زمانی بررسی تغییرات پدیده‌های متفاوت. وجود الگوریتم‌های پیش‌فرض, دسترسی به دو نوع داده‌های ple‏ همراه با 
تصحیح اتمسفری و همچنین داده‌های مربوط به محصولات و پروداکت‌های مختلف. عدم نیاز به سخت‌افزارهای 
قدرتمند و کاهش حجم داده‌ها اشاره نمود. مطالعات بسیاری در زمینه تحلیل و بررسی حشکسالی با استفاده از 


داده‌های سنجش دور صورت گرفته است که به‌عنوان نمونه در ادامه به تعدادی از آن‌ها اشاره می‌گردد. 


۲ - مواد و روش‌ها 
۱-۲ - محدوده موردمطالعه 
محدوده موردمطالعه در اين پژوهش شهرستان داراب یکی از شهرستان‌های استان فارس می‌باشد که در جنوب 
شرق اين استان قرار دارد و بين ۵۶ درجه و ۰٩‏ دقيقه و ۰۵ ثانيه تا ۵۵ درجه و ۲۷ دقيقه و ۰۷ asl‏ طول شرقی تا ۲۸ 
درجه و ۰۱ دقيقه و ۲۷ ثانيه تا ۲۸ درجه و OT‏ دقيقه و ۳۳ ثانيه عرض شمالی و به فاصله ۲۵۰ کیلومتری از شهر 
The Standardized Precipitation Index (SPI)‏ 1 


2 Terra 
3 Google Earth Engine (GEE) 


شيراز به‌عنوان مركز استان فارس قرار دارد (شكل ۱). به لحاظ اقليمى اكثر بارش‌های شهرستان داراب در طى دوره 
سرد سال روى مىدهد و در تابستان‌ها به دليل دورى از ale‏ رطوبتى و نبودن عامل صعود هواى gh ys‏ حالت 
هواى گرم و خشك كويرى يبدا می‌کند. بر اساس داده‌های ثبت شده توسط ايستكاه سينويتيك حسن‌آباد موجود در 
شهرستان ميانكين بارش سالانه آن طى دوره آماری ۲۰ ساله (۱۳۹۹ — ۱۳۷۹ ۲۳۷۳ ميلىمتر و ميانكين دماى هواى 
ثبت شده آن ۲۲/۲ درجه سلسیوس می‌باشد. شهرستان داراب از نظر ژئومورفولوژی جزء زاگرس جنوب‌شرقی 
محسوب می‌شسود که در اين منطقه» زاگرس به حداکثر پهنای خود رسیده و چهره ظاهری ناهمواری‌ها با ساير 
بخش‌های زاگرس متفاوت است (علائی طالقانی. ۱۳۸۲). میانگین ارتفاع شهرستان داراب حدود ۱۶۹۱ متر از سطح 
دریاهای آزاد می‌باشد. قسمت‌های زیادی از شهرستان به علت وجود سنگ چینه‌های شور کننده مانند سنگ‌های 
سازند هرمز سازند آغاجاری و سازند ساچون از لحاظ بهره‌برداری با محدودیت مواجه می‌باشند. شهرستان مذکور از 
لحاظ پوشش گیاهی بسیار غنی و دارای انواع مختلف گیاهان خودرو می‌باشد و تنها در برخى از قسمت‌ها که از 
خاک‌های شور تشکیل شده گیاهان شوری پسند مانند كز دیده می‌شود (سلیمانی‌شیری» ۱۳۹۱). 
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شكل -١‏ تصوير ماهوارهاى موقعيت جغرافيايى شهرستان داراب 


FF 


54°26'1S"E 


-Y-Y‏ داده‌های اقليمى و ماهواره‌ای مورداستفاده در پژوهش 

جهت انجام يزوهش حاضر از دو نوع osla‏ استفاده شده است. یکی داده‌های ميزان بارش روزانه ثبت شده توسط 
ايستكاه سینوپتیک حسن آباد داراب» به‌عنوان تنها ايستكاه سينويتيك موجود در شهرستان كه داراى آمار طولانى بارش 
ثبت شده می‌باشد و داده‌های آن از سال ۱۳۷۸ تا سال ۱۳۹۹ جهت محاسبه فاش SPI‏ در بازه‌های زمائی مختلف 


مورداستفاده قرار گرفت. آمار مربوطه از سازمان هواشناسی کشور دریافت شد و ضمن انجام آزمون‌های آماری مانند 


۱۸۰ 


جغرافیا و مخاطرات محیطی 


شمارة چهلم 


آزمون همگنی و آزمون داده‌های پرت. فرایند بازسازی مقادیر ناقص و مفقودی از روش نسبت‌ها بر روی آن‌ها انجام 
شد. داده دیگر مورداستفاده در اين مطالعه شامل تولیدات ماهواره‌ای پوشش گیاهی (NDVI)‏ و دمای سطح زمین 
(LST)‏ مربوط به سنجنده مودیس ماهواره تررا می‌باشد که فرایندهای ارزیابی و استفاده از آن‌ها توسط سامانه تحت 
وب گوگل ارث انجین صورت گرفت. شکل (۲) تصویری از محيط کدنویسی سامانه گوگل ارث انجین را نشان 


می‌دهد. با توجه به اينكه پدیده خشکسالی به‌صورت آنی رخ نمی‌دهد و تأثیرات OF‏ بعد از گذشت چندین ماه و يا 


چندین سال نمایان می‌شود بنابراین در اين پژوهش به‌منظور ایجاد سری زمانی تغییرات از تصاوير مربوط به ماه‌های 


خرداد تا شهریور هرسالء به Wo‏ اينكه در اين زمان پوشش گیاهی در وضعیت مناسبی قرار دارد. در طی دوره آماری 


۰ ساله (۱۳۹۹ — ۱۳۷۹) و درمجموع ET‏ تصوير با تو جه به همخوانی زمانی و همچنین تبدیل تصاویر ۸ روزه به 
۱۹ روزه» مورد تجزیه‌وتحلیل جهت انجام پایش و ارزیابی خش كسالى مورداستفاده قرار گرفت. در ادامه جدول A)‏ 


ویژگی‌های تصاویر مورداستفاده شده را نشان می‌دهد. 


6 * ه 9 © هه و و Ki‏ + 


on 


> Imports (1 entry) Q 


Search places and datasets... 


© د + 


@ code.earthengine.google.com 


Google Earth Engine 


[scrips Docs Assets 0۳ I مه ] مدوم‎ _-[] ۰-- -1۲ ۷۰ 1621 inspector (oy Mase 


a Use print(...) to write to this console. 


۱ var table: Table users/hesmaeili191/darab 


var ndvi = ee. ImageCollection("MOOIS/@06/M001301" ) 


CE: | 


Filter scripts. 


k Drought_TCI_VCI gao 


~ Owner (1) 1 Map.centerobject(table); 
cay TS Rit 2 nap د ا‎ (tobie); » ImageCollection HODIS/006/H0D13Q1 (8 elenents) ason | 
Sade 3° /* Object (@ properties) 
4 IIIIII year ++ vel 2500 
i Drought & LST modis 5 // modis vegetation product 
6 
7 
8 


ISON 


* Image (1 band) 


-filterBounds(table) 
-filterDate(' 2000-05-21", '2000-09-21") 


&LST_MODIS_day_night 


lisST_chloraNOAAGFSOP25 || و‎ :select( ۱0( Object (@ properties) 
۱ test! 107 .map(function(img){ Mean_VCI JSON 
k ECMWF_ERAS_LAND_HOURLY 11 ver cier هت‎ 
12 return clip.multiply (0.6001) 
ee 1 13 -copyProperties(img,['system:time_start’, "system:time_end']); ۱ ImageCollection MODIS/@@6/MOD11A2 (15 elements) ISON 
sea_surface_elevation_temp_salt 24 زر‎ 
k sentinel 15 X 
k sse_salt 1 | 6 ۱ ImageCollection MODIS/@06/MOD11A2 (15 elements) ISON v 
نيا‎ a p ۳۳ Pt 
0 9 2 ۲ 0 Layers a jpu 
soe 
+ ¥ 850 ١ فیروزآناد فرانیسد‎ j rt 
يوشكات‎ ee خبره‎ ۷ La 
بات‎ 
|  ومروخ‎ 80 جو اسعیدفه‎ 
| aes 
۳ اسن ون قير ددم‎ Myre 
vgs 0 دوي أن حاحي ا قير‎ m 
3 / 
sin “ws Em a 
رخاشاب‎ ry 
m © ٠ deb كمي‎ 
م‎ 53 Loy 
500016 SS تت 20218 ندر‎ jas Wiss دنه‎ S/S 


شكل ۲ - محيط كدنويسى سامانه كوكل ارث انجين (GEE)‏ 


حدول -١‏ مشخصات تصاوير استفاده شده در طى سال‌های ۱۳۷۹ تا ۱۳۹۹ 


قدرت تفکیک 


نام ماهواره نام سنحنده 


نام محصول 


كد محصول 


تررا موديس دماى سطح زمين MOD11A2 LST)‏ ۰ متر ۸روزه 
تررا موديس يوشش كياهى MOD13Q1 (NDVI)‏ ۰ متر ٩‏ روزه 


1 Land Surface Temperature (LST) 


سال دهم بررسی سری زمانی خشکسالی کشاورزی شهرستان داراب .... ۱۸۱ 


۳-۲ - شاخص‌های استفاده شده حهت بررسی خشکسالی 
۱-۳-۲ - شاخص وضعیت دمایی ‘TCD‏ 
این شاخص جهت تعیین دمای مرتبط با تنش پوشش گیاهی که ممکن است ناشی از رطوبت بیش‌ازحد نیز باشد 
بە‌صورت رابطه (۱) تعریف می‌شود TOTS)‏ 1940 وست و همکاران ؛ ۲۰۱۹): 
Ty. -T‏ 


TCI =| عل‎ 


Pua I mn ale) 


كه در رابطه فوق T‏ دمای سطح زي Tua Tun LST)‏ به ترتيب كمينه و بيشينه دماى سطح زمين در 
محدوده موردمطالعه می‌باشد. 

۲-۳-۲ - شاخص وضعیت پوشش گیاهی VCD‏ 

شاخص وضعيت يوشش كياهى از طريق رابطه (Y)‏ زیر محاسبه می‌شود EÉ S)‏ 440 وست و COL Raw‏ 
0۳۰۹ 

_{ NDVI-NDVI,, 

NDVI „ -NDVI „, 
يفيه‎ Rls Su NDV nn 5 57 اد قار ت رک كام‎ a 5 کان‎ 


رابطه (۲) 


و كمينه شاخص NDVI‏ می‌باشد. اين شاخص شرایط غالب پوشش گیاهی را نشان می‌دهد. به‌طوری‌که اگر مقدار آن 
به صفر نزدیک شود خحشكسالى شدید رخ می‌دهد (فرامرزی و نوری. ۱۳۹۶: (VEO‏ مطابق با مطالعات و بررسی‌های 
انجام شده مقادیر شاخص‌های TCL VCI‏ در ينج كلاس به شرح جدول (Y)‏ قابل دسته‌بندی می‌باشند (رستمی و 
همکاران, ATA‏ 114 سلیمانی و همکاران. ۱۳۹۸: ۲۱؛ اسدی می‌آبادی» ۱۳۹۶: Te‏ رضایی‌بنفشه ۱۳۹۶: ۱۱۸؛ ساطانی‌فر و 


همكاران. ۱۳۹۸: NEV‏ تتکبیل و گامگی" ۲۰۰6 فرهمند " ۲۰۱۹): 


جدول ۲- کلاس‌های مختلف شاخص‌های TCI VCI‏ 


وضعیت خشکی ارزش پیکسل وضعیت BES‏ ارزش پیکسل 
خشكى بسيار شديد 1< خشکی ملایم ۰ ۳۰ 
خشكى شدید 2-۰ ۱۰ بدون خشكى >t‏ 


1 Temperature Condition Index (TCI) 
2 Kogan 

3 West et al 

4 Kogan 

5 West et al 

6 Thenkabail & Gamage 

7 Farahmand 


\AY‏ جغرافیا و مخاطرات محيطى شمارة چهلم 


وضعيت خشكى ارزش پیکسل وضعيت خشكى ارزش پیکسل 
خشكى متوسط Ye ye‏ = = 


۳-۳-۲۳ - شاخص بارش استاندارد شده SPI‏ 

جهت برآورد خشکسالی هواشناسی» شاخص SPI‏ که توسط مک‌کی و همکاران (۱۹۹۳) ارائه شده و دارای 
قابليت و اعتبار بالايى Lak e‏ مورداستفاده بسيارى از محققين و پژوهشگران قرار می كيرد. شاخص SPI‏ در بازه 
زمانى بلندمدت با وقوع سيلابهاء تغييرات سطح آب موجود در مخازن سدها و همجنين نوسانات منابع آب 
زیرزمینی در ارتباط می‌باشد (مک‌کی و همکاران KORGAN e‏ جهت محاسبه شاخص SPI‏ از رابطه )۳( استفاده مى شود: 


sp )۳( رابطه‎ 
SD 


اق او Malay‏ زتره لد عه ail‏ ار مات بارش ا by‏ مرس اش ایا 
ميانكين بارش درازمدت ايستكاه موردنظر و SD‏ انحراف معيار بارش می‌باشد. طبقات اين شاخص به‌صورت 
جدول (۳) می‌باشد. مقادير درجهبندى اين شاخص بدین‌صورت است كه مقادير مثبت Sly‏ مقادير بیشتر از ميانه 
SS‏ و اوور بور رانک کش پاش لقانت لت جنا عم كلق hyo so‏ ا sie‏ 


اښست: 


جدول ۳- طبقهبندى خشكسالى با استفاده از شاخص SPI‏ 


ترسالى بسيار شديد 
ترسالى شديد 
ترسالى متوسط 
ترسالى خفيف (نرمال) 
خشكى خفيف (نرمال) 
خشكى متوسط 
خشكسالى شديد 


خشكسالى بسيار شديد 


SPI مقدار‎ 
۲ > SPI 
۱/۵ < SPI 1/44 
۱/۹ > SPI < \ 
۰/۹۹ > SPI > ۰ 


۰/۰۰ > SPI < -۵ 


-۱/۰۰ > SPI > -1/24 


-۱/۵ > SPI < -1/44 


SPI < ۰ 


1 McKee et al 


سال دهم بررسی سری زمانی خشكسالى کشاورزی شهرستان داراب .... ۱۸۳ 


۳ - نتایج و بحث 

VCI و‎ TCI بررسی وضعيت خشكسالى شهرستان داراب بر اساس شاخص هاى خشكسالى‎ -۱ -Y 

ازجمله شاخصهاى ماهواره‌ای به كار رفته در اين پژوهش جهت پایش و ارزيابى خشكسالى شاخص وضعيت 
دمايى (TCD‏ می‌باشد که جهت محاسبه آن از يروداكت دماى سطح زمين LST)‏ استفاده شد. Y pal‏ از شاخص TCI‏ 
جهت تعيين دماى مرتبط با تنش يوشش گیاهی و همجنين تنش ناشى از رطوبت بیش ازحد استفاده می‌شود. اگر 
مقدار عددى ناشى از اين شاخص به ۰ برسد نشان‌دهنده عدم خشكسالى می‌باشد. از طرف ديكر هرجه مقدار اين 
شاخص به صفر نزدیک شود نشان‌دهنده شرايط خشکی و وقوع خشکسالی می‌باشد؛ بنابراين با توجه به ميزان 
مساحت طبقات مختلف اين شاخص در اين پژوهش مشخحص شد که بالاترین مساحت خشکی با طبقه خیلی شدید 
در سال ۱۳۷۹ به میزان ۲۳۵/۶۳ کیلومترمربع و پس‌ازآن به ترتیب در سال‌های ۱۳۸۹ ۱۳۹۷ ۱۳۹۰ و ۱۳۹۵ بیشترین 
مساحت خشکی خیلی شدید در شهرستان داراب اتفاق افتاده است. از سویی دیگر کمترین میزان مساحت خشکی 
سال‌هایی که بیشترین مساحت بدون خشکی را داشته‌اند مربوط به سال‌های ۱۳۹۸ و ۱۳۹۹ است که هر دو سال 
مذکور به‌صورت مشابه و یکسان دارای مساحتی برابر با ۶۲۶/۶۱ کیلومترمربع می‌باشند و در مرتبه بعد سال ۱۳۸ با 
مساحت ۲ کیلومترمربع بیشترین مساحت بدون خشکی را داشته است. مساحت pls‏ طبقات شاخص TCI‏ 
در جدول )£( به تفکیک سال ارائه شده است. در ادامه نقشه‌های حاصل از شاخص TCT‏ برای دوره موردبررسی به 
تفکیک هرسال در اشکال (۳ و )٤‏ ارائه شده است. از یافته‌های قابل‌توجه پژوهش اين است که نقشه‌های شاخص 
خشكسالى TC]‏ برای دو سال ۱۳۹۸ و ۱۳۹۹ شرایط کاملاً یکسان و مشابهی را نشان می‌دهند به‌طوری‌که قبلاً نيز 
بیان شد مساحت طبقات مختلف خشكى در آن‌ها كاملا یکسان می‌باشد که اين امر حاکی از اين است که شرایط 


تنش دمایی و همچنین تنش ناشی از رطوبت در طی اين دو سال در محدوده موردمطالعه یکسان بوده است. 


شکل ۳- نقشه‌های شاخص وضعیت دمایی (TCI)‏ و وضعیت پوشش گیاهی (VCI)‏ شهرستان داراب در سال 
۱۳۹۹ 


۸٤‏ جغرافيا و مخاطرات محيطى شمارة چهلم 


BS 
Ss 


“Ovo 25 50 Ok” a = — im wel: 
5 EAHA بدون خشکی خشکی ملایم خشکی متوسط خشکی شدید خشکی خيلى شديد‎ 


شكل -é‏ نقشههاى شاخص وضعيت دمايى (TCI)‏ شهرستان داراب برای سال‌های ۹ نا ۱۳۹۸ 


حدول ۶- مساحت طبقات خشكسالى شهرستان داراب در شاخص وضعيت دمایی (TCD‏ به كيلومترمربع 


سال خشکی بسیار شدید خشکی شدید خشکی متوسط خشکی ملایم فاقد خشکی 
۱۳/6۹ ا ۸۰۹ ۹/1۱ 0747 Yoww‏ 
VAAY \ eY o/y Waa 10۹/۰۲ ۱۳۸۰‏ ۳۵ ۱۳۸ 
Yovv/14 ۹/1 ۱۱۸ ۱/۳۹ \AVEN ۱۳۸۱‏ 
VEA MANY ۸۱/۰ 1/10 ۱۳۸۲‏ 01 
TITTY ۱۳/۳۳ 47/7۹ Wer) AE ۱۳۸۳‏ 
YVAL/W Ves VAV ۹77 AY£/\ \ \EL/AV ۱۳۸۶‏ 
FAVA ۱۳/۸ Ve YEA OVO (2 ۱۳۸۵‏ 
۱۳۸۹ ام( YAVA/AY A۸1 0/771 TAY‏ 
YALV/V4 1/71 7/۰۸ Vey ۱۳۹/۲ AY‏ 
RAA ۹۹/4 ۱۰۸۹ 1/۶ \A\/00 NYM‏ 


سال دهم 


سال 

۱۳۸۹ 
۱۳۹۰ 
۱۳۹۱ 
۱۳۹۲ 
۱۳۹۳ 
۱۳۹۶ 
۱۳۹۵ 
۱۳۹۹ 
۱۳۹۷ 
۱۳۹۸ 


۱۳۹ 


خشکی بسیار شدید 
۳۳/۸ 
YAY‏ 
\0V/ Y‏ 
AA 0A‏ 
۹/۳ 
WALE‏ 
ANE‏ 
VEI‏ 
۲10/۰۱ 
۱۱۱/۳۷ 


(۵2 


بررسی سری زمانی خشکسالی کشاورزی شهرستان داراب .. 


خشکی شدید 
VEE‏ 
۰2۳/۹ 
۱۱۱۱۸۳۹ 
Vor‏ 
0۹0/7 
E/E‏ 
111/07 
۵( 
WAY‏ 
LAY‏ 


WAY 


خشكى متوسط 
۹۸0/0۹4 
۱/۷ 
۹4۸/4۷ 
11/1۰ 
|0/۷ 
tvy‏ 1 
10/۲ 
9/۸۷۶ 
۹۳۸/۳۷ 
۸۳۹۲ 


۸۳۹۲ 


٩/۲ 
۱۳ 
۱۰۳۹/۹۹ 
10/0۰ 
۱۰۳۹/۰۸ 
ANA 
27/11 
40710 
MWA 


MIA 


فاقد خشكى 
۱۳۹ 
۱۳۸۹۸ 
۱۳۷۵۷۰۱ 
۱۳۸۱۷۸ 
YMI‏ 
TVAN/ET‏ 
۸1/41 
ی 
۱۳۹۷۷۶ 
AA)‏ 
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با توجه به شاخص وضعيت پوشش كياهى (VCD‏ و بررسىهاى صورت گرفته بيشترين مساحت خشكى با 
طيقه خشکی خيلى شديد مربوط به سال ۱۳۹۲ با مساحت ۲۸۰/۸۰ کیلومترمربع و در مراتب بعدى سال‌های ۳۸۳ 
۷۲ و ۱۳۹۱ بيشترين مساحت خشکی با طبقه خيلى شديد را تجربه كردهاند. كمترين ميزان مساحت طبقه خشکی 
خيلى شديد نيز مربوط به سال ۱۳۹۹ به مساحت ۱۱/۹۹ كيلومترمربع می‌باشد. از سويى دیگر سال 1145 با مساحت 


۸ کیلومترمربع بيشترين مساحت طبقه بدون خشكى را دارا بوده است. در مراتب بعدى سالهاى ۱۳۹۹ ۱۳۹۸ 


و ۱۳۹۷ بيشترين مساحت طبقه بدون خشکی را بر اساس شاخص VCI‏ تجربه كردهاند كه اين امر بيانكر وقوع 


سال آرائه شده اسث. 


جدول ۵- مساحت طبقات خشكسالى شهرستان داراب در شاخص وضعيت يوشش كياهى (VCI)‏ به 
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۱۸۹ جغرافيا و مخاطرات محيطى شمارة جهلم 


سال خشکی بسیار شدید خشکی شدید خشکی متوسط خشکی ملایم فاقد خشکی 
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شكل ۵- نقشه‌های شاخص وضعيت پوشش گیاهی (VCD)‏ شهرستان داراب برای سالهاى ۱۳۷۹ تا ۱۳۹۸ 


سال دهم بررسى سرى زمانى خشكسالى كشاورزى شهرستان داراب .... Wwy‏ 


۲-۳ - بررسى وضعيت خشکسالی با استفاده از شاخص SPL‏ در ايستكاه سينويتيك محدوده مورد مطالعه 

همانگونه كه در قسمت‌های قبل عنوان شد در اين پژوهش از داده‌های ايستكاه سینوپتیک حسن آباد داراب كه 
داراى آمار طولانى مدت بارش ثبت شده می‌باشد. از سال ۱۳۷۹ تا سال ۱۳۹۹ جهت محاسبه شاخص خشكسالى 
SPI‏ در بازه‌های زمانى يك ماهه. سه ماهه» شش ماهه و دوازده ماهه استفاده شد. جهت محاسبه شاخص ٩۳]‏ از 
نرم‌افزار MATLAB‏ استفاده شد. بر اساس يافتههاى تحقيق بيشترين ميزان خشکسالی بر اساس شاخص مذكورء در 
سال ۱۳۷۹ همراه با طبقه خشکسالی خفيف و سال ۱۳۸۷ همراه با طبقه خشكسالى متوسط روى داده است. در ادامه 
جهت تعيين مطلوبترين شاخص ماهوارهاى نشان دهنده خشكسالى در منطقه مورد مطالعه به مقايسه بين مقدار 
عددى شاخص SPI‏ برای هر يك از سال‌های موردبررسی با ميزان مساحت‌های فاقد خشكسالى سال متناظر حاصله 
از دو شاخص TCI‏ و VCI‏ به دست آمده از تصاوير ماهوارهاى اقدام شد شكل CV)‏ (شايان ذكر است جهت يكسان 
سازى واحد مساحت مناطق فاقد خشكسالى حاصل از شاخص TCI‏ و VCI‏ با شاخص SPT‏ از روش استاندارد 
سازى Z-Score‏ استفاده شده است). همجنين مطابق با بررسی صورت كرفته بيشترين مقدار ضريب همبستگی 
پیرسون به ميزان ٠/76‏ بين SPI‏ دوازده aal‏ و شاخص VCI‏ مشاهده كرديد كه اين مطلب نمایانگر شاخص 


ماهواره‌ای بهینه و مطلوب نشان‌دهنده وضعیت حشکسالی در منطقه مورد مطالعه می‌باشد. 


مقابسه شاخص SPI‏ در بازه زمانی ۱۲ و ۶ ماهه با مساحت ola‏ بدون خشکسالی تصاویر ماهواره 
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شكل *- نمودار مقايسه شاخص SPI‏ با مساح تهاى بدون خشكى حاصل از شاخص‌های ماهوارهاى TCI‏ و 
VCI‏ 


ملا جغرافیا و مخاطرات محيطى شمارة چهلم 


٤‏ — جمع بندى 

در طى سالهاى اخير پیشرفت‌های صورت گرفته در زمينه سنجثرازدور و تصاوير ماهواره‌ای موجب شده تا 
داده‌های بسيار زيادى با قدرت تفکیک مکانی و زمانى مناسب و همجنين قابليت پوشش مناطق وسيع در دسترس 
پژوهشگران محققين و همجنين دست‌اندرکاران علوم محيطى و علوم زمين قرار كيرد. ازجمله زمينههاى كاربردى 
استفاده از اين داده‌های ارزشمندء ارزيابى و يايش خث كسالى مناطق مختلف به‌ویژه مناطقی می‌باشد كه در آن‌ها به 
دليل شرايط خاص خود امكان برداشت و ثبت داده‌های حاصل از مشاهده و اندازهكيرى زمینی آن‌ها وجود نداشته و 
يا اينكه بسيار سخت» صعبالعبور و يرهزينه است كه شهرستان داراب از اين امر مستثنى نیست؛ لذا در اين پژوهش 
جهت ارزيابى و Gul‏ سرى زمانى خشكسالى كشاورزى شهرستان داراب طى يك دوره ۲۰ ساله (۱۳۹۹ — ۱۳۷۹) 
از داده‌های ماهوارهاى و همجنين شاخصهاى پرکاربرد سنجش‌ازدوری از قبيل شاخص‌های وضعيت دمايى (TCT)‏ 
و شاخص وضعيت يوشش گیاهی (VCD‏ استفاده شد. نتايج ارزيابى شاخص TCI‏ حاكى از اين است كه در طى 
دوره موردبررسى بيشترين مساحت خشكى با طبقه خيلى شديد در سال ۱۳۷۹ به ميزان ۲۲۵/۶0 کیلومترمربع و پس 
از آن به ترتیب در سال‌های ۸۱۳۸۹ ۰۱۳۹۷ ۱۳۹۰ و ۱۳۹۵ يشترين مساحت خشكى طبقه خيلى شديد در شهرستان 
داراب اتفاق افتاده است. همچنین نيز سال‌هایی که بیشترین مساحت بدون خشکی را داشته‌اند مربوط به سال‌های 
۸ و ۱۳۹۹ که هر دو سال مذکور بصورت مشابه و یکسان دارای مساحتی برابر با ۶۲۱۶/۶۱ کیلومترمربع می‌باشند 
و همچنین نیز در مرتبه بعد سال ۱۳۸۲ با مساحت ۲۳ کبلومترمربم وسعت بیشترین ميزان طبقه بدون خشكى 
را داشته است. یکی از یافته‌های قابل توجه پژوهش حاضر اين است که نقشه‌های شاحص خشکسالی TCI‏ برای دو 
سال ۱۳۹۸ و ۱۳۹۹ شرایط کاملاً یکسان و مشابهی را نشان می‌دهند كه اين امر حاکی از اين است که شرایط تنش 
دمایی و همچنین تنش ناشی از رطوبت در طی اين دو سال در محدوده مورد مطالعه یکسان بوده است. 

با توجه به شاخص وضعیت پوشش گیاهی (VCD‏ بیشترین مساحت خشکی با طیقه خشکی be‏ شدید 
مربوط به سال ۱۳۹۲ با مساحت ۲۸۰/۸۰ کیلومترمربع و در مراتب بعدی سال‌های ۰۱۳۸۳ ۱۳۸۲ و ۱۳٩۱‏ بیشترین 
مساحت خشکی با طبقه Ube‏ شدید را تجربه کرده‌اند. کمترین میزان مساحت طبقه خشکی خیلی شدید بر اساس 
شاخحص VCE‏ نيز مربوط به سال ۱۳۹۹ به مساحت ۱۱/۹۹ کیلومترمربع می‌باشد. از سویی دیگر سال ۱۳۹۹ با 
مساحت ۱۱۳/۷۸ کیلومترمربع بیشترین مساحت طبقه بدون خشکی را دارا بوده است. در مراتب بعدی سال‌های 
۵۹ ۱۳۹۸ و ۱۳۹۷ بیشترین مساحت طبقه بدون خشکی را بر اساس شاخص VCI‏ تجربه کرده‌اند. جهت تعیین 
مطلوب‌ترین شاخحص ماهواره‌ای ULE‏ دهنده خشكسالى در شهرستان داراب به مقایسه بين مقدار عددی شاخص 
SPI‏ برای هر یک از سال‌های دوره موردبررسی با میزان مساحت‌های فاقد خش كسالى سال متناظر حاصله از دو 
شاخص TCI‏ و VCI‏ به دست آمده از تصاوير ماهوارهاى اقدام شد. مطابق با بررسی صورت گرفته بیشترین مقدار 


ضريب همبستكى پیرسون به ميزان ۲ بین SPI‏ دوازده ماهه و شاخحص VCI‏ مشاهده گردید که اين مطلب 


نمایانگر شاخص ماهوارهاى بهينه نشان‌دهنده وضعيت خش كسالى در منطقه مورد مطالعه مىباشد و اين امر داراى 
همسويى و مشابهت با پژوهش‌های انجام شده توسط سالازار' و همكاران (۲۰۰۸» رضايى بنفشه و همكاران (۱۳۹۶) 
و اسدی می‌آبادی و همكاران (۱۳۹۶) در محدوده‌های موردبررسى خود می‌باشد كه در مطالعه خود جهت بررسى 
وضعیت خشکسالی از تصاویر ماهواره‌ای و شاخحص VCI‏ استفاده نموده‌اند و به اين نتيجه رسیدند كه شاحص VCI‏ 
دارای قابلیت زیادی جهت پایش خش كس الى کشاورزی مناطق و همچنین در رابطه با عملکرد محصولات زراعی 
می‌باشد. یکی از مزیت‌ها و تفاوت‌های مهم اين مطالعه نسبت به تحقیقات مشابه نامبرده و ساير تحقیقاتی از اين قبيل 
که با کمک تصاویر ماهواره‌ای به بررسی مخاطره خشکسالی مناطق مختلف اقدام نموده‌اند. اين است که در اين 
مطالعه از تعداد زیادی تصاوير ماهواره‌ای در قالب سری زمانی )£19 تصوير به فاصله هر ١6‏ روز) جهت بررسی 
E ELET‏ استفاده شسده اس که اسفاده از فرایك سری زمانی تصاویرو افزایش تعناد: تضاوترن باك 
می‌شود كه يايش و بررسی تغییرات خشکسالی در سال‌های مختلف بهتر نمود بيدا AS‏ بنابراین اين امر موجب شد 
که تحقیق حاضر از اين حيث نسبت به مطالعات نامبرده و ساير مطالعاتی که تنها یک تصوير منتخب از یک سال را 
به‌عنوان نماینده و الگوی سال مورد نظر لحاظ می‌کنند. دارای برتری باشد. مزایای دیگر تحقیق صورت گرفته نسبت 
به ساير مطالعات مشابه اين است که تحقیق حاضر در محبط سامانه تحت وب گوگل ارث انجین و به‌وسیله فرایند 
کدنویسی به‌عنوان فرایندی پویا و با قابلیت‌های بالا انجام شده است که کاهش بسیار زياد زمان فرایندهای پردازش 
تصاوین امکان ايجاد سری زمانی بررسی تغییرات پدیده‌های متفاوت. دسترسی به دو نوع داده‌های خام و همراه با 
تصحیح اتمسفری و همچنین داده‌های مربوط به محصولات و پروداکت‌های ماهواره‌ها و سنجنده‌های مختلف. عدم 
نياز به سيستم سخت افزاری و نرم‌افزاری قدرتمند و همچنین نیز کاهش حجم داده‌ها از قابلیت‌های مهم و اساسی 
سامانه گوگل ارث انجین و فرایند کدنویسی آن محسوب می‌شود که استفاده از اين مهم برتری تحقیق حاضر نسبت 
به ساير تحقیقات مشابه را نشان می‌دهد. از دیگر نتایج برجسته تحقیق حاضر اين است که استفاده از داده‌های 
سنجش‌ازدور و سامانه گوگل ارث انجین جهت پایش و بررسی خحشکسالی در مناطقی که فاقد داده‌های مشاهداتی 
حاصل از برداشت‌های زمینی می‌باشد بسیار سودمند می‌باشد لذا به مدیران و برنامه‌ریزان حیطه علوم محیطی 
شهرستان داراب توصیه می‌شود که در اهداف و اقدامات آتی خود جهت برنامه‌ریزی مواردی مانند مقابله با مخاطره 
خشکسالی به اين مهم توجه نمایند. 


1 Salazar 
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